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INTRODUCTION 


In order to predict the convective heat- transfer distribution for 
the windward surface of the Space Shuttle entry configuration, one must 
develop engineering correlations for the three-dimensional, compressible 
boundary- layer. Since the aerodynamic heating rates generated by a tur- 
bulent boundary-layer may be several times greater than those for a laminar 
boundary-layer at the same flight condition, the correlations must include 
a transition criteria suitable for the complex flow fields. Because the 
windward surface of the Orbiter is composed of a large number of thermal 
prot-'ction tiles, the transition criteria must include the effect of the 
distributed roughness arising from the joints and possible tile misalignment. 

During tests in which a ring of spherical roughJiess elements were 
located in a supersonic flow past a cone. Van Driest and Blumer (ref. 1) 
observed variations in the relative roles played by the disturbances in the 
basic flow field and those resulting from the p»*esence of roughness elements. 
For some conditions, the disturbances associated with the basic flow field 
were predominant in establishing transition, whereas for other flows, the 
roughness elements dominated the transition process. 

However, the correlation for the effect of roughness is complicated 
when other transition-related parameters interact. Morisette (ref. 2) 
found that although the effective roughness Reynolds number increases sig- 
nificantly in the presence of a favorable pressure gradient near the center- 
line, much smaller roughness was required to promote transition near the 
'houlder of an Orbiter configuration, where again there was a favorable 
pressure gradient (this one associated with cross flow). McCauley et al 
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(ref. 3) found that the spherical roughness elements required to trip the 
boundary layer on sphere noses were several times larger than the boundary- 
layer thickness, whereas the trips required for a cone were within the 
boundary layer. Heat- transfer data (ref. 4) obtained in Tunnel B of the 
Arnold Engineering Development Center (AEDC) for an 0.04-scale Orbiter 
indicated that a ring of spherical trips, which were 0.079 cm. (0.031 in.) 
in diameter and were 0.11L from the nose, caused the transition location to 
move considerably upstream of the natural transition location (i.e., that 
for a smooth body). In the same test program (ref. 4), a simulated inter- 
face gap between two insulation materials, which was 0.102 cm. (0.040 in.) 
wide by 0.203 cm. (0.080 in.) deep and was located at x = 0.02L, had no 
measurable effect on boundary-layer transition at o = 40® and Re , = 

*,L 

8.6 X 10®. In a series of tests using delta-wing Orbiter models (ref. 5), 
premature boundary-layer transition was observed on a model having siri..ilated 
heat-shield panels with raised joints. Slot joints, however, did not cause 
premature transition of the boundary layer. The former model featured a 
series of transverse panels 0.635 cm. (0.250 in.) wide separated by a raised 
retaining strip 0.025 cm. (0.010 in.) wide by 0.0025 cm. (0.001 in.) high. 

The panels on the model with slotted joints were 0.635 cm. (0.250 in.) square 
separated by slots 0.020 cm. (0.008 in.) wide by 0.005 cm. (0.002 in.) deep. 
The Reynolds number (Re^ 1^) for these tests ranged from 6.5 x 10® to 
9.0 X 10® using a model 0.403 m. (1.321 ft.) long. 

The stability of laminar boundary-layers has been found to be 
significantly affected by heating or by cooling (usually indicated para- 
metrically by a temperature ratio, or enthalpy ratio, such as T /T or 
T^/Tg). Lees (ref. 6) found that heat-transfer from the fluid to the wall 
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stabilized a laminar boundary layer for two-dimensional disturbances and 
that. If there Is a sufficient amount of cooling, the boundary layer on 
a hydrodynamically, smooth configuration could be completely stabilized. 
Reshotko (ref. 7) presents data supporting the trend toward complete 
stabilization. Using the flight data of Rumsey and Lee (ref. 8), Reshotko 
observed that, for conditions outside the predicted region of complete 
stabilization, transition did occur (as expected), but at relatively 
high Reynolds numbers. 

However, the unqualified prediction that cooling stabilizes the 
boundary layer cannot be made since transition reversal has been observed 
by numerous workers as the model Is "cooled" (e.g., refs. 9 and 10). As 
a result of cooling the boundary layer, there Is a relative Increase In 
the magnitude of the disturbance due to a fixed roughness (ref. 11). A 
further complication Is associated with wind tunnel data. Since a low 
value of Tyij. may be obtained either by cooling the wall or by heating 
the test gas, alternative effects may arise. Using data from a single 
tunnel, Wagner et al (re;. 12) noted that reducing by heating the 

flow significantly decreased Re ^ , possibly because of nonunlform mix- 
ing of the supply gas in the stagnation chamber. 

Data from an experimental program which was conductef^ to Investigate 
what effect tile misalignment representative of a reasonable manufacturing 
tolerance has on heat transfer and transition criteria In the plane-of- 
symmetry of the Shuttle Orbiter have been analyzed (e.g., ref. 13). The 
vertical tile misalignment simulated on the 0.0175-scale model was approx- 
imately 0.1451 cm. (0.0571 in.) full-scale. Furthermore, the surface 
temperature for the Tunnel B tests was essentially constant at 0.42T^. As 
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noted in ref. 13, the presence of tile misalignment did not significantly 
affect the transition locations over the range of test conditions considered. 
Another test program was conducted in Tunnel F where the windward surface was 
roughened by a grit blasting technique. The surface temperature for the 
Tunnel F tests varied from 0.14T^ to 0.2bT^. At the higher Reynolds numbers 
of the Tunnel F tests, the transition location moved near the nose. The 
roughness elements became large relative to the boundary layer and became 
effective as tripping elements. However, there were not sufficient parametric 
variations to establish suitable correlations. 

The present report presents calculations which indicate the effect that 
variations in the flow-field model would have on correlation of wind-tunnel 
data. The theoretical heat-transfer distributions are compared with experi- 
mental heat-transfer distributions obtained in Tunnel B at the Arnold Engineer- 
ing Development Center (AEDC) using a 0.0175 scale model of the Space Shuttle 
Orbiter Configuration for which the first 80% of the windward surface was 
roughened by a simulated tile misalignment. The experimental heat-transfer 
data were used to determine the transition locations. Data were obtained for 
a Mach number of 8 over a Reynolds number range (based on model length) from 
1.862 X 10® to 7.091 x 10® with surface temperatures from 0.114T^ to 0.435T^. 
with tile-misalignment heights of 0.0025 cm. (0.0010 in.) and of 0.0051 cm. 
(0.0020 in.). 
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NOMENCLATURE 

Cp pressure coefficient (eqn 1); or specific heat for a constant 

pressure process (eqn 2) 

h metric scale-factor used to represent the radius of an equi- 

valent body-of-revolution (see Fig. 7); or the local heat 
transfer coefficient, q/(T -T ) 

I W 

H stagnation enthalpy 

ht heat-transfer coefficient for the stagnation point of the 

reference sphere 

k height of the misaligned tiles 

L axial model length, 0.5734 m. (1.881 ft) 

M Mach number 

p local stat’c pressure 

q local heat- transfer rate 

R gas constant 

^ref radius of the reference sp^^ere, 0.00533 m. (0.0175 ft.) 

RCns Reynolds number based on flow conditions behind a normal shock 
(eqn 6) 

Re^ Reynolds number based on local flow properties Integrated along 
the wetted distance along a streamline (eqn A) 

RCg Reynolds number based on local flow properties and the momentum 
thickness 

Re^^l^ Reynolds number based on free-stream flow properties and the 
model length (eqn 5) 

s entropy (eqn 2); or wetted distance along a streamline (eqn 4) 
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T temperature 

^init temperature used in correlation of transition location (see 
eqn 14b) 

recovery temperature 

U streamwise component of the velocity 

X axial coordinate 

y coordinate measured normal to the model surface 

a angle of attack 

•k 

6 displacement thickness 

6 momentum thickness 

y viscosity 

p density 

relative transition location defined in eqn 17 

?2 relative transition location defined in eqn 18 

Subscripts 

e evaluated at the edge of the boundary layer 

i evaluated at the particular run of interest (see eqn 17) 

ns evaluated downstream of a normal shock wave 

tr evaluated at the transition location 

w evaluated at the wall 

» evaluated at the free-stream conditions 



EXPERIMENTAL PROGRAM 




A primary objective of the present Investigation was to determine what 
effect tile misalignment representative of a reasonable manufacturing toler- 
ance has on the heat-transfer distribution In the plane of symmetry of the 
Shuttle Orbiter. To do this, the heat- transfer data from three test pro- 
grams conducted in Tunnel B of the AEDC were studied. The data were obtained 
at a free-stream Mach number of 8 over a range of free-stream Reynolds number 
(based on model length) from 1.862 x 10® to 7.091 x 10®. The surface tem- 
perature was varied from 0.1 14T^ to 0.435T^. The data presented in the pre- 
sent report were obtained at angles-of-attack of 30® and 40®. The test con- 
ditions for these data are summarized in Table 1. 

Models 

The basic model used in the test programs (see the sketch of Fig. 1) was 
a 0.0175-scale model of the Space Shuttle configuration defined by Rockwell 
drawing VL70-000139 and designated ^todel 29-0. Twenty seven coaxial surface 
thermocouples were used to obtain the heat-transfer-rate distribution for the 
windward plane- of- symmetry. The locations of the heat-transfer gages in tfe 
windward pi ane-of- symmetry are indicated in Fig. 1. To study the effect of 
tile misalignment, selected tiles were precisely etched (or deposited, as 
discussed in ref. 14, depending on the misalignment height) on the windward 
surface, so that they were slightly above the model surface. The misaligned 
tilfs formed a herringbone pattern (symmetric about the plane of symmetry) 
covering the windward surface of the Orbiter model up to the tangent line 
of the chines from x = 0.02L to 0.80L. The raised tiles, which were selected 
randomly, represented 25 % of the tilsn in the area of Interest, as shown in 
the photograph of Fig. 2. The selectf.<; tiles were 0.267 cm. (0.105 in.) square. 
The model surface for each of the three test programs is summarized below. 
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0H4A test program . - For the 0H4A tests, the surface of the 0.0175-scale 
Orbiter model was smooth. Thus, the transition locations determined from 
the heat- transfer distributions serve as the reference, or smooth-body, 
transition locations. The reader Is referred to re*. 15 for the basic 
data and for additional Information regarding these tests. 

MH2A test program . - For the MH2A tests, the misaligned tiles were deposited 
to a height of approximately 0.0025 cm. (0.0010 in.). The vertical misalign- 
ment, thus simulated, was 0.1451 cm. (0.0571 in.) full-scale. Thus, the 

-5 

nondimensionalized tile height (k/L) for this model was 4.43 x 10 . The 
basic data for this tile misalignment height (designated k, in the present 
report) are presented in ref. 16. An analysis of the heat- transfer data and 
the transition locations was presented in ref. 13. 

MH2B test program . - For the MH2B tests, the surface surrounding the tiles 
was removed until the misaligned tiles were approximately 0.0051 cm. (0.0020 
in.) in height. This misalignment, which is designated k^ in the present 
report, corresponds to a full-scale vertical misalignment of 0.2903 cm. 
(0.1143 in.). For the kg misalignment, k was 8.86 x 10~^L. Additional 
information about the model and the basic data for the MH2B test program is 
presented in ref. 17. 



THEORETICAL FLOW FIELDS 


Theoretical solutions for a nonsimilar, laminar boundary- layer were 
generated for the pitch plane of the Orbiter model at angles-of-attack of 
30° and of 40°. These theoretical boundary-lay'' ' solution^ were computed 
using a modified version of the code described in ref. 18. Required as 
input for the code are the flow conditions at the edge of the boundary 
layer, the radius of the "equivalent" body-of-revolution, and the wall- 
temperature distribution. The metric scale-factor describing the stream- 
line divergence was used to represent the radius of the equivalent body-of- 
revolution in the axisymmetric analogue for a three-dimensional boundary- 
layer. Two different flow-field models were used to generate the required 
inviscid solution for the plane of symmetry. The pressure distribution and 
the streamline-divergence characteristics for the first flow model, desig- 
nated "Mod-Newt NSE", represent modified Newtonian flow. For this flow 
model, it was assumed that the fluid at the edge of the boundary layer was 
that which hac assed through the normal portion of the bow shock wave and 
had accelerated isentropically from the stagnation point to the local 
static pressure which was defined by: 

Cp = Cp^ sins ( 1 ) 

The streamwise distributions of the static pressure, of the entropy at the 
edge of the boundary- layer, and of the cross-sectional radius for the 
equivalent body-of-revolution for the second flow model, designated "Var 
Ent", were computed using the code described in ref. 19. Because the bow 
shock-wave is curved, the entropy varies throughout the shock layer (the 
magnitude of the entropy change across the shock wave depending upon the 
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free-strean Nach nunber)- The entropy at the edge of the boundary layer 
depends not only on the streaiMise location but on the local boundary-layer 
thickness and» hence, on the Reynolds nunber. However* the variation of 
the local entropy at the edge of the boundary layer due to the variations 
In the boundary-layer thickness was believed to be of second order Importance 
for the present study. Therefore* the local entropy was assumed to be that 
obtained In the theoretical solution for Re , = 1.9 x 10® at the Orbiter 
surface. 

Surface static-pressure distributions for the plane of synmetry are 
presented in Fig. 3. For o * 30®* the theoretical distributions for the 
two flow models are compared with the experimental pressures presented In 
ref. 20. The theoretical values for the "Var Ent" flow model are In very 
good agreement with the experimental values. The theoretical pressures 
for the "Mod-Newt HSE" are significantly below the experimental values. 

Once the local static pressure and the local entropy have been defined, 
the local Mach number at the edge of the boundary layer can be calculated 
using the following relations: 



Note, that In all calculations, the air has been assumed to behave as a 
perfect gas. The Mach number distributions presented in Fig. 4 clearly in- 
dicate that the "Var Ent" flow model provides a more realistic representation 
of the actual flow field than does the "Mod-Newt NSE" model Although exper- 
imental values were not available for a = 40®, extrapolation of the data for 
o = 30® and for a = 35® Indicates that the "Var Ent" flow model should pro- 
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vide a reasonable representation of the actual flow field at an angle-of- 
attack of 40“ . 

The Hach number distribution across the laminar boundary- layer is pre- 
sented in Fig. 5. The theoretical distribution obtained using the Univer- 
sity's NONSIMBL code is in excellent agreement with that obtained using the 
BLIMP code as calculated by Goodrich et al (e.g., refer to ref. 19). Ex- 
perimental values of the Mach nun^er which were obtained at similiar test 
conditions (ref. 21) are included for comparison. The experimental values 
are in fair agreement with the two sets of theoretical calculations. 

Streamwise distributions of the theoretical value of the Reynolds number 
at the edge of the boundary layer are presented in Fig. 6. The local values 
of Re^, where 

s 

r p u 

*^0 ' 

have been divided by the free-stream Reynolds number based on model length. 

The local Reynolds numbers were calculated using the local entropy distribu- 
tion for the "inviscid" flow-field for Re . = 1.9 x 10®. Since the entropy 

»,L 

is evaluated at the edge of the local boundary layer, which is a function of 
the free-stream conditions, the value of the local entropy would also depend 
on Re I . If we had accounted for the entropy variation normal to the surface, 
the nondimensionalized Reynolds number at the edge of the boundary layer would 
have depended on the free-stream Reynolds number. This Reynolds-nunber ef- 
fect was neglected for the purposes of the present report. Note that the 
local Reynolds number which was calculated using the Var Ent flow model is 
roughly twice the corresponding value for the Mod-Newt NSE flow model (refer 
to Fig. 6a). Thus, both the local Mach number (refer to Fig. 4) and the 
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local Reynolds number are significantly greater for the Var Ent flow field 
model. 

The metric which represent the cross-sectional radius of the equivalent 
body-of-revolution is presented as a function of the axial coordinate in Fig. 

7. The metric coefficients were calculated using the relations described 
by De Jarnette (ref. 22) and by Rakich and Mateer (ref. 23). The streamlines 
near the plane of symmetry diverge much more rapidly for the Var Ent flow- 
model. 

With the exception of the Mach number profiles for a laminar boundary- 
layer which were presented in Fig. 5, the theoretical values presented thus 
far describe the inviscid flow field (specifically, conditions at the edge 
of the boundary layer). Since the streamwise distributions for the entropy 
at the edge of the boundary layer and for the static pressure were provided 
by Or. W. D. Goodrich of the Johnson Space Center (JSC), NASA, one would 
expect that the values of parameters such as the edge Mach number generated 
using the NONSIMBL code would correlate well with the values calculated by 
others. Thus, to substantiate correlations based on the viscous portion of 
the code, the heat-transfer-rate distribution calculated using the University's 
NOMSIMBL code for one of the test conditions was compared with the distribution 
calculated at the JSC using the BLIMP code. As shown in Fig. 8, the rates 
calculated using the NONSIMBL code are within 15%, or less, of the values 
calculated using the BLIMP code. 

The effect of surface temperature on the heat-transfer rates in the plane 
of symnetry is illustrated by the calculations presented in Fig. 9. The 
results are presented as the local heat-transfer coefficient divided by the 
heat-transfer coefficient for the stagnation point of a 0.00533 m. (0.0175 
ft.) radius sphere as calculated using the theory of Fay and Riddell (ref. 
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24). For the present report, the recovery factor r has been set equal to 
unity. Distributions of the dimensionless heat-transfer-coefficient ratio, 
h/ht j are presented for each of the flow models for the extreme values 
of wall temperature covered in the experimental program. For each case, 
the surface temperature was assumed to be constant along the entire length 
of the model. For each of the three flow models, the nondimensionalized 
heating was slightly greater for the colder surface. This should be ex- 
pected since the boundary layer is thinner at the lower surface temperature, 
the velocity gradients and, hence, the energy generated within the boundary 
layer by viscous dissipation are greater. 

Because of the streamwise variation in the heat-transfer, the surface 
temperature did not remain constant during an experimental run as was as- 
sumed in the calculations presented in Fig. 9. Therefore, the theoretical 
heat-transfer distribution which was calculated assuming the wall tempera- 
ture to be constant (T^ = 0.423 T^) was compared in Fig. 10 with that cal- 
culated using the experimentally-determined surface temperatures (the sur- 
face temperatures varied from T^ = 0.450 T^ to T^ = 0.429 T^). These the- 
oretical heat-transfer distributions are compared with the corresponding 
experimental heat-transfer coefficients. However, the improvement is 
relatively small. Therefore, subsequent solutions for the laminar boundary- 
layer assume the surface temperature to be constant along the entire length 
of the model. The values used in computing theoretical solutions will 
represent "nominal" values from the experimental program. 



EXPERIMENTAL DATA 
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The heat- transfer measurements are presented as a dimensionless ratio 
of '.eat-transfer coefficients, h/h^ which involves the ratio of the 
rieasured, local heat-transfer rate to the theoretical heat-transfer rate to 
the stagnation point of a 0.00533 m. (0.0175 ft.) radius sphere as calculated 
usitg the theory of Fay and Riodell (ref. 24). For purposes of data pre- 
sen; ation, the recovery factor has been set equal to unity. Other para- 
meters used in the data correlations include the free-stream Reynolds number 
based on model length. Re , , where 

“>,L 

p U L 

^ <R> 00 • _ . 


a id L is the model length, 0.5734 m. (1.881 ft.). The Reynolds number behind 

a nofmel shock ^ is 
ns 


Re 


ns 


= ^ns ^ns ^ref 
^ns 


where r^^^ = 0.00533 m. (0.0175 ft.). 


( 6 ) 


Heat-Transfer Distributions for a = 30° 

The eirect of the height of the misaligned tiles on the heat-transfer 
distribution for the plane-of-symmetry when the Orbiter is at an a of 30° 
i illustrated in the data presented in Fig. 11. Data are presented from 
tests where the surface temperature of the model was approximately 300°K 
(540°RV «:or the geometries tested in the wind tunnel and for this wall 
teniD .rature, tile misalignment did not significantly affect the heat-transfer 
rates in regions where the boundary layer was either laminar or turbulent. 
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Furthermore, the boundary layer remained laminar along the entire length 
of the model for the lower Reynolds number, i.e.. Re , = 1.92 x 10®. At 

“’,L 

the higher Reynolds number, i.e., Re^ = 7.03 x 10®, the heat-transfer 

measurements for the model with k * 4.43 x 10*®L departed from the laminar 

distribution at approximately the same location as observed for the smooth 

-5 

body. Increasing the tile misalignment height to k = 8.86 x 10 L caused 

the transition location to move slightly upstream and caused a more rapid 

streamwise increase in heating, i.e., the heating increased very rapidly 

with distance in the transitional flow regime. 

Theoretical heat- transfer distributions for a laminar boundary layer 

when = 0.416 T,. and Re , = 3.67 x 10® are included in Fig. 11 for com- 
w t «,L ^ 

parison with the experimental distribution. Since (as has been noted 
already) the variations in the local entropy at the edge of the boundary 
layer due to the variations in the boundary- layer thickness had been neglect- 
ed, the theoretical value for the dimensionless ratio of the local heat- 
transfer coefficient was independent of the test Reynolds number. Solutions 
are presented for both flow models, i.e., the "Var Ent" and the "Mod-Newt 
NSE" flow models. For x > 0.2L, the Var Ent model provides a better predic- 
tion of the laminar measurements. It is interesting to note that whereas the 
theoretical solution describes the laminar boundary layer adjacent to a 
smooth surface, it provides a reasonable prediction of the heating rates to 
the tile-roughened surface. 

The data presented in Fig. 11 are from tests for which the model surface 
was at room temperature. Thus, the wall temperature was approximately 42 
percent of the stagnation temperature of the tunnel flow. This is signif- 
icantly higher than the values for the ratio T /T. which would be enco ntered 

W 

in flight. Therefore, to more nearly simulate the flight values, the model 
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was cooled by exhausting liquid nitrogen over the model prior to injecting 
it into the free stream. Cooling the model (and, hence, the boundary layer) 
causes the boundary* layer thickness to decrease. Theoretical Mach-number 
profiles are presented for the /arious surface temperatures of the test 
program in Fig. 12 for Re^ = 1.9 x 10^ and in Fig. 14 for Re^ = 7.1 x 10^. 
Shadowgraphs of the corresponding flow fields are presented in Figs. 13 and 
15. At the lower Reynolds number, the theoretical solutions (which neglect 
any perturbation of the flow due to the presence of the misaligned tiles) 
indicate that the tile is submerged within the subsonic portion of the 
boundary layer. Although weak waves generated by the tiles appear in each 
shadowgraph (see Fig. 13), they only become easily visible for the coldest 
surface temperature, i.e., = 0.128 T^, for which the boundary layer is 

thinnest. Although the presence of the perturbation waves indicates at least 
a weak flow-field disturbance, the heat-transfer rates in the laminar region 
do not appear to be significantly affected by the tile misalignment (see 
Fig. 16a). However, whereas the boundary layer remained laminar along the 
entire length of the model for T = 0.423 T. and for T = 0.340 T^, the 
heat-transfer distribution indicates that the boundary layer becomes transi- 
tional for x > 0.70L when was decreased to 0.128 T^. Thus, as a result 
of cooling the boundary layer, there is a relative increase in the magnitude 
of the disturbance due to a fixed roughness. The effect is obviously a non- 
linear function of surface temperature. 

The Mach number profiles at the higher Reynolds number, i.e.. 

Re I = 7.1 x 10®, indicate that the misaligned tiles would protude into 
the sonic regions of an unperturbed boundary layer at x = 0.255L (see Fig. 

14). Thus, one would expect that the presence of the misaligned tiles would 
alter the basic flow field. That the misaligned tiles significantly perturb 
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the flow field is verified in the shadowgraphs of Fig. 15. Easily visible 
waves emanate from the tiles at all surface temperatures. At the lower 
surface temperatures, the highly vortical flow which results when the turbu- 
lent boundary layer bounces over the misaligned tiles can be seen in the 
photcg'-aphs, see Fig. 15(c) and 15(d). As a result, transition occurs well 
upstream cn the fuselage, as evident in the heat- transfer distributions of 
Fig. (16c). For = 0.231 T^, the heating increases from x = 0.08L to 
X = O.IOL, and then decreases nearly following the laminar theory for *-he 
gage at x = 0.1 5L. Downstream of x = 0.1 5L, the heating increases rapidly 
with distance as the boundary layer becomes fully turbulent. The presence 
of a strongly favorable pressure gradient may have inhibited the transition 
process in the region O.IOL < x < 0.15L. For T = 0.171 T., the heat-trans- 

W U 

fer measurements depart from the laminar distribution at x = 0.07L. However, 

there is little difference in the transition location for T = 0.401 T* and 

w t 

that for = 0.313 T^. These data indicate that the primary impact of 
cooling was to shorten the transition length for these two test conditions. 
Thus, again any correlation between the transition location and the surface 
temperature would be nonlinear. 

The effect of surface temperature on the heat-transfer distribution is 
illustrated in the data of Fig. 16. The theoretical solutions indicate 
th c thinning the boundary layer by surface cooling increased the non- 
dimensional ized value of the local heat-transfer coefficient. The heating 
increase is attributed to the increased gradients and viscous dissipation 
which result when the boundary layer thickness is decreased. At a given 
location, the nondimensional ized heat-transfer coefficient measured for the 
coldest surface temperature was more than for the highest surface temperature, 
as oredicted bv theory. However, the values obtained at the intermediate 
surface temperatures did not fall between these extremes. Thus, although 
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there are considerable differences in the experimental values in the laminar 
region for the various surface temperatures, there is no consistent trend 
with surface temperature. However, the coaxial surface thermocouple gages 
used to measure the heating rates were not calibrated to show variations of 
this magnitude. 


Heat-Transfer Distribution for a = 40° 

The effect of the height of the misaligned tiles on the heat- transfer 
distibution for the pi ane-of- symmetry when the Orbiter is at an a of 40° 
is illustrated in the data presented in Fig. 17. The measurements are taken 
from tests where the surface temperature of the model was approximately 
300° K (540°R). For the geometries tested in the wind tun 1 and for this 
wall temperature, tile misalignment did not significantly affect the heat- 
transfer rates in regions where the boundary layer was either laminar or 
turbulent. Increasing the tile misalignment height to k = 8.86 x 10 L 
caused the transition location to move only slightly upstream and caused 
a more rapid streamwise increase in heating, i.e., the heating increased 
very rapidly with distance in the transitional flow regime. 

The decrease in the boundary- layer thickness due to cooling the model 
is illustrated by the theoretical Mach number profiles presented in Fig. 

18. At this Reynolds number, the theoretical solutions (which neglect 
any perturbation of the flow due to the presence of the misaligned tiles) 
indicate that the tile is submerged within the subsonic portion of the 
boundary layer. However, as the wall temperature is decreased, the velocity 
profile is such that the flow is transonic near the top of the tile. 
Snadowgraphs of the corresponding flow fields are presented in Fig. 19. 
Although weak waves generated by the tiles appear in each shadowgr'ph, they 
only become easily visible for the coldest surface temperature, i.e., T = 
0.138 T^, for which the boundary layer is thinnest. 
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The effect of surface temperature on the heat- transfer distribution is 
illustrated by the data presented in Fig. 20. As was the case for the 
a = 30° measurements, these heat- transfer gages did not verify the theoretical 
trend with surface temperature. There is a definite effect on the transition 
location. For both values of the free-stream Reynolds number, the transi- 
tion location moves rapidly upstream as the surface is cooled, since there 
is an increase in the relative magnitude of the disturbance. 


Transition Correlations 

The heat-transfer distributions were used to determine the point at 
which boundary-layer transition occurred in the plane of symmetry. The 
experimentally-determined transition location was that "point" at which 
the heat transfer first deviated from the laminar distribution. The transi- 
tion locations, thus determined, are presented in Fig. 21 and in Table 1. 

The transition location is a function of the angle of attack, the free- 
stream conditions (specifically, the Reynolds number, since the Mach number 
is essentially constant for these Tunnel B tests), the surface temperature, 
and the misalignment height. Thus, for a given angle of attack, the transi- 
tion location is a function of the dimensionless variables: k/L, and 




ns* L’ 


(7) 


Although the use of a local Reynolds number would be desirable, Re^^ was 
chosen because it was both simple to calculate and yet could account for 
the variations in the free-stream Mach number. Since data were obtained (1) 
."or the smooth body over a range of Reynolds number, (2) over a range of 
tile height for a single surface temperature, and (3) over a range of sur- 
face temperature with a fixed tile height, it was possible to develop an 
empirical correlation for the transition location as a product: 
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I I = (Smooth-Body Correlation) (Roughness 

corr (8) 

Factor) (Temperature Factor) 

This empirical correlation is developed in the following paragraphs using 
the data for an angle-of-attack of 30®. 

Smooth-body correlation . - Referrina to the experimental Iv-determined 
transition locations for the smooth-bodv Orbiter model, i.e.. k^. with 
T^ = 0.42 T^ at an angle-of-attack of 30®, it is clear that: 

Smooth-Body Correlation = ^ ®1 ^°^10 ^®ns 

A least-squares fit of the data yields the numerical values for a^ and 
a^: 

Smooth-Body Correlation = 4.3131 - 0.9858 109io 

Equation (10) provides a good correlation of the data, as can be seen in 
Fig. 22. 

Roughness factor . - As was discussed in reference 13 and as is illustrated 

in the data of Fig. 21, the transition locations for a tile-roughened 

-5 

model where k = 4.43 x 10 L were essentially the same as those recorded 

for a smooth body. However, when the heights of the misaligned tiles 
-5 

were 8.86 x 10 L, the transition locations moved forward approximately 
0.05 L, which is the distance between adjacent heat-transfer gages. 

Using another set of data obtained in Tunnel B, Hube (ref. 25) observed 

that, "Projecting tiles were effective boundary- layer trips". The height 

-5 

k of the transition-promoting tiles varied from 9.69 x 10 L to 

-5 

48.45 X 10 L. This corresponds to a full-scale tile misalignment of 
0.3175 cm. (0.1250 in.) to 1.5875 cm. (0.6250 in.). The effect of 
tile-misalignment height in promoting transition is clearly nonlinear. 
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The roughness factor which represents the forward shift in transi- 
tion due to the misaligned tiles was assumed to be given by: 

. 2 

Roughness Factor = 1.0 - b^^. (^) (11) 

The boundary-layer perturbation due to tiles of a given height should 
increase as the boundary layer thins, i.e., as the Reynolds nuiti)er in- 
creases. To represent this Reynolds-number dependei.ee, b^^ was assumed 
to be a linear function of Re^^. A least-squares fit of the experiment- 
ally-determined transition locations for T^ = 0.42 T^ and for the k^, , 

and surface finishes (which are given as the first eighteen test con- 
ditions of Table 1 for an angle-of-attack of 30°) was used to establish 
a value of b^^ for each of the Reynolds numbers tested. Another least- 
squares fit was used to correlate these values of b^^.. The resultant 
roughness-factor correlation is: 

fi t 2 

Roughness Factor = 1.0 - [9.3805 x 10° + 174.94 Re^^] (-^) (12) 

The transition locations aiven bv equation (12) are comoared with the 
exoerimental locations for T = 0.42 T^. in Fia. 23. The emoirical cor- 
relation is in aood aoreement with the individual exoerimental values. 

Temoerat ' ire factor . - Data were obtained for the k^ surface finish, i.e.. 
k = 8.86 x 10 ^ L, over a range of surface temperature from 0.11 T^ to 
0.42 T^. These data were used to determine the temperature factor, 
which represents the forward shift in transition due to cooling the 
boundary layer. Thus, the temf 'rature factor represents the ratio of 
the transition location observed at the temperature of interest divided 
by the transition location observed at the reference temperature, which 
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is 0.42 for the present tests. The temperature factor: 




k,^,0.42T^,Rej^g 


is presented in Fig. 24 as a function of the surface temperature. Cool- 
inq the wall to 0.32 did not significantly affect the transition loca- 
tions. However, when the wall was cooled to 0.22 or to 0.12 
matic forward movement of the transition location was observed. Fur 
more, the effect was very sensitive to Reynolds number. For the lower 
surface- temperatures, i.e., = 0.12 or 0.22 T^, transition occurred 

very near the nose (depending on the Reynolds number), approaching a 
minimum transition length. Thus, the expression for the temperature 
factor should be such that: 

(1) for large surface temperatures, i.e., as T^ approaches 0.42 T^, 
the temperature factor approaches unity, 

(2) for the lowest surface temperatures , the temperature factor ap- 
proaches a minimum value, and 

(3) the temperature factor changes rapidly with surface temperature 
once the temperature is below a critical minimum. 

Data provided by Dr. Goodrich (ref. 26) indicated that cooling the wall 
from 0.42 T^ to 0.12 T^ had little effect on the transition location for 
the smooth body and had only a moderate effect on the transition location 
for the model with the surface finish. The general form of the pro- 
posed temperature factor is: . /T - T 

^2 I w init 

"mr\ ’’t 

Temperature Factor = 1.0 - C-j e ' 
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where Is the wall temperature at which the transition location 

reaches its upstream-most location. Decreasing the surface temperature 
below this value would not move transition forward. Because the tile- 
induced boundary- layer perturbation is a function of the height of the 
tile relative to the displacement thickness, the temperature-dependent 
parameter should reflect an interdependence on the Reynolds number. Thus, 
^init ® function of the Reynolds number, A correlation of the experi- 
mentally-determined transition locations yields the following expression 
for the temperature factor: 

3.322 X 10‘^ r .v ■ ^init\ 

rk7LT~ “t; — 

Temperature Factor = 1.0 - 0.85 e ' ' (14a) 

where 

c 2 

T^-nit " ^ 

(see Fig. 24). The empirical correlation for the temperature factor is 
compared with the experimental data in Fig. 24. Although the least- 
squares expression for fits the individual calculations quite well, 

there are significant di-^ferences between the correlation values and the 
individual experimental values for the temperature factor in the range 
0.11 0*23 T^. These differences will be reflected in the com- 

parison between the experimental and the correlation values for the 
transition location. 

Resultan t correlation for jUie a - 30° . Tunnel B data . - When the three 
factors, as given by equations (10), (12), and (14), are substituted into 
equation (8), the empirical correlation for the t ansition location for 
the tile-roughened Orbiter model at an angle of attack is 30® is: 



10 ■'^ns 



where T- - = 82.033 + 4.0573 x lO”® (Re„)^. The transition location 
imc 

moved significantly forward for a highly cooled surface, for high 
Reynolds nunbers, and for relatively large tile-misalignment height. 

These parameters indicate that premature transition occurs when the 
height of the misaligned tiles was of the order of the displacement 
thickness. Note that the correlation uses only parameters which repre- 
sent the test conditions, i.e., •'/L. It does not re- 

quire knowledge of the local flow properties and, therefore, does not 
depend on the flow model used to calculate the local flow field. 

The transition locations calculated using this empirical correlation 
are compared with the experimentally-determined transition locations in 
Table 1. For extreme test conditions, two limiting considerations were 
employed. For the highest surface temperature and for the k^ and the k^ 
surface finishes, the correlation yields transition locations in excess 
of 1.00 at the lowest Reynolds number. Values in excess of unity, i.e., 
x^^ > L, were interpreted as representing wholly laminar boundary-layers. 
For that case where T^ was less than (refer to Fig. 24 and note 

the hexogonal symbol), the temperatur actor was assumed to be equal to 

its lower limit, 0.15. This avoided the error which would result if the 
actual temperatures were used, since the exponent would then be positive 
and the correlation of Fig. 24 violated. 



In addition, the difference between the en^irical value and the 
experimental value of the transition location for a given test condition 
is presented in Fig. 25 as a function of Re^^. For the majority of the 
test conditions, the two values are within 0.05 L of each other. Since 
the heat- transfer gages were 0.05 L apart over most of the fuselage (see 
Fig. 1), the difference is approximately equal to the ability to define 
the transition location. The greatest differences result in those tests 
where the surface had been cooled considerably. For these tests, the 
transition location was very sensitive to 

reflected the sensitivity of the displacement thickness relative to the 
tile height on these parameters. As noted earlier, even though the 
least-squares expression for fit the individual calculations quite 

well, there were significant differences in the values for the temperature 
factor at these conditions. Thus, the correlation is considered very good 
for these Tunnel B data. 

Applicability of the correlation to other data. - When a correlation to 
predict the transition location for a particular configuration is developed 
using data from a single wind tunnel, one must ask how general is such a 
correlation? Thus, the correlation of the Tunnel B data, represented by 
equation (15), was compared with transition locations determined using 
heat-transfer distributions obtained on a similar model in Tunnel F (see 
refs. 27 and 28). por the Tunnel F tests, the Mach number varied from 
10.73 to 12.06 while the free-stream Reynolds number (based on model 
lenqth) varied from 1.17 x 10® to 17.63 x 10®. The surface temperature 
varied from 0.16 T^ to 0.27 T^. The surface was roughened over the 
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80t of the windward surface using a grit-blasting technique. The average 
peak-to-valley distance for ten readings in a 0.25 on. (0.10 in.) length 
as read from a photo-micrograph was 0.0041 cm. (0.0016 in.). As can be 
seen in the data oresented in Fio. 26. surface roughness did not have a 
significant effect until the Reynolds nimber was sufficiently high to 
Cause transition near the nose. At these high Reynolds numbers, the 
roughness elements became large relative to the boundary layer and be- 
came effective as tripping elements. 

As can be seen in Fig. 26, the experimentally-determined transition 
locations for the smooth-body tested in Tunnel F fall significantly below 
the smooth-body correlation developed using the Tunnel 8 data, i.e., 
equation (10). However, by changing the value of the constant a^, one 
obtains an approximate fit of the smooth-body Tunnel F data. Specifical- 
ly, the relation 



4.175 - 0.9858 log^Q Re^^^ 


(161 


provides a reasonable correlation of the Tunnel F data. The "prematute'’ 
transition of the boundary layer which occurred in the Tunnel F tests may 
have been due to tunnel noise, to nonuniform mixing of the relatively 
high- temperature flow, or to some other factor. It is also possible 
that a correlation of the transition location as a linear function of the 
log^Q Re^g is not appropriate for these data and that this parameter 
does not satisfactorily relate data from two different tunnels. 

Thus, although the transition locations calculated using equation 
(15) are in good agreement with the individual experimental values from 



Tunnel 6, they do not agree with the Tunnel F data. Although the various 
factors, which are given in equations (10), (12), and (14), qualitatively 
describe the effect of the various parameters, caution should be used 
before applying equation (15) to predict the transition location for a 
Shuttle Orbiter exposed to significantly different flow conditions, i.e., 
those of another wind tunnel or of flight. 

Correlations Based on Local Flow Properties 

The Var Ent flow model was used to calculate the flow properties at 
the edge of the boundary layer and the momentum thickness at the 
experimentally-observed transition location. The ratio (R®0/Mg)^^» 
thus calculated, is presented in Fig. 27 as a function of the nondimen- 
sionalized transition location, x^^/L. For those runs with the higher 
surface temperatures, i.e., 0.32 T^ or 0.42 T^, the transition locations 
were not "»uch different from the reference- temperature, smooth-body 
transition locations. For these runs 



both for an angle-of-attack of 30° and of 40°. As often noted previously, 
cooling the surface increased the tile-induced flow perturbations and 
caused transition to move forward. As a result, the value of 
decreased continuously, approaching 100. 

The previous discussion of the experimentally-determined transition 
locations has underlined the fact that the tile-induced flow perturbations 
become strongest when the height of the misaligned tiles is of the order 
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of the displacement thickness. The transition location moved signifi- 
cantly forward as surface cooling at the higher Reynolds numbers thinned 
the boundary layer. The transition location for the perturbed flows 
relative to the corresponding locations on the reference configuration 
are presented as a function of ii. Figs. 28 through 30 and of 

5*/k in Flos. 31 through 33. Since smooth-body data were obtained only 
for ungl'is-of-attack of 25**, 30®, and 35* (ref. 15), the reference 
transition-locations for an a of 40® were those for the k.j- tile-roi jhened 
model (ref. 16). Recall that there was no consistently measurable differ- 
ence between the transition locations for the smooth-body (k^) and those 
for the k.j- tile roughened model (k * 4.43 x 10”^ L) for an a of 30® (see 
Fig. 21). Thus, the relative transition locations ^or o = 30® are pre- 
sented as the ratio where 


(v\ 
(— ) 


07) 


For a = 40®, the relative transition locations re presented as the ratio 
where 


^2 



•^i * ^wi ’*^®ns 


,0.42T^,Re|^^ 


(18) 


These nondimensionalized relative transition-locations represent the ratio 
of the experimentally-determined transition location at the tile misalign- 
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ment height of a particular run (k.)» at the surface temperature of that 
run at the Reynolds number of that run divided by the transi- 

tion location measured at the same Reynolds number but at the reference 
surface finish (either or k^) and at the reference surface temperature 
(0.42 T^). 

The values of 6 / Re^/k and 6*/k at x = 0.1 L and x = 0.2 L were 
calculated for a nonsimilar, laminar boundary layer. For some conditions, 
boundary- layer transition occurred upstream of these x-locations. Thus, 
the actual value of & for these cases would be significantly different 
than the calculated value. Theoretical solutions were generated for the 
Mod Newt NSE flow model at a = 30®, for the Var Ent flow model at a = 30® 
and for the Var Ent flow model at o = 40®. As has been discussed, the 
theoretical flow models do not reflect tile-induced perturbations. All 
the test conditions of Table 1 are not represented in Figs. 28 through 
33. Theoretical solutions were not generated for the Mod Newt NSE flow 
model for the a = 40® flow conditions or for the k^ -surface finish at an 
a of 30®. Furthermore, at the lowest surface temperature, i.e., T = 

0.12 T^, the displacement thickness calculated at x = 0.1 L using the 
Mod Newt NSE flow model was often negative or very small (depending upon 
the Reynolds number). Thus, the relative transition location is repre- 
sented by the range of values for these runs and is presented at an 
approximate value of «*/k. For a given surface temperature, the para- 
meter 6*/ Re^/k is essentially constant over the range of Reynolds number. 
Thus, for that surface temperature and a given surface roughness, the 
data tend to cluster in a group when this parameter is used as the 
abcissa (see Figs. 28 through 30). 



For either nondimensionalized 6*-parameter, the relative transition- 
location z is near unity (i.e.> the tile-induced perturbations have a 
negligible effect on the transition location) at the highest values of 
the 6*-parameter. As the 6*-parameter is decreased (i.e., the boundary 
layer is cooled, the test Reynolds number is increased, or the tile- 
misalignment height is increased), transition moves gradually forward, i.e., 
£ decreases slowly. In this region, i.e., 0.8 < c < 1.0, the magnitude 
of z does not exhibit any dependence on the Reynolds number except through 
the value of the 6*-parameter. Below a critical value of the 6*-parameter, 
the relative transition location decreases rapidly (i.e,, transition moves 
rapidly upstream toward the nose) as the boundary layer thins. Below the 
critical value of the 6*-parameter, the Reynolds number has a strong effect 
on the magnitude of c. Note that the critical value of the fi*-parameter 
evaluated at a given location is independent of the flow-field tnocel used 
in the calculation of 5* and Re^, and is independent of the angle-of-attack 
for these wind tunnel results . 

The relative transition locations were correlated as a function of 
the theoretical, laminir values of 6*/k at x * O.IL using a linear, 
least-squares fit of those data for which 6* > 0.75 k. The least-squares 
correlations are essentially identical for the o = 30® data. I.e., for 
the Mod Newt NSE model: 


^1 


0.8404 + 0.0489 



O.IL 


(19) 


while for the Var Ent model: 

= 0.8397 + 0.0517 



( 20 ) 


The least-squares correlation for the a « 40® results is: 
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C, = 0.9017 + 0-0074 (211 

* O.IL 

Recall that the heat-transfer distributions and, therefore, the transi- 
tion locations were not measured for the smooth-bodv (k^l at a » 40®. 

It was assumed that the transition locations for the k-j -tile-roughened 
model would have been the same as those for the smooth bodv. Thus, they 
Were used as the reference transition-location in evaluating the £2 
parameter. As a result, there were no measurements where £ was unity, 
as had been the case for the o = 30® correlations (see Figs. 31 and 32). 

The absence of such points is the reason why eqn. (21) is significantly 
different than eqn (19) or eqn (20). Nevertheless, eqn (19) or eqn (20) 
would provide a reasonable correlation of the £2 nieasurements for 6* > 

0.75 which are presented in Fig. 33. 

The relative transition locations for both angles-of-attack are pre- 
sented in Fig. 34a as a function of the theoretical, laminar values of 
6*/k at X = O.IL for both flow models. Note that for 6* ^ 0.75k, the ex- 
perimental values of £ >^0.81, i.e., the roughness-perturbed transition 
locations are within of the reference, smooth-body transition loca- 
tions. At the lowest temperatures, 6* was so thin, (i.e., 6* < 0.75k) 
that the presence of the misaligned tiles moved transition well upstream. 
Thus, 0.75k is the critical value of 6* at x * O.IL. Using a linear, 
least-squares fit of those data for which 6* at x * O.IL was greater than 
0.75k yields the correlation: 

£ = 0.8489 + 0.0468 ( f-\ (22) 

V / X = O.IL 

The parameter £ is used, since the correlation fits both the £-|- and the 
£ 2 -values presented in Figs. 31-33 and, therefore, is valid for both angles 
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of-attack and for both flow models. Using this correlation, when the 
theoretical value of the displacement thickness at x = O.IL is equal to 
0.75k, i = 0.884. Thus, the “nominal" transition location moves upstream 
by just over 10% when 6* at this station is of the order of the tile mis- 
alignment height. 

The relative transition locations are presented in Fig. 34b as a 
function of the theoretical, laminar values of 5*/k at x - 0.2L for both 
flow models. Note that for 6* ^ 1.50k, the experimental values of e 1 
0.81. The linear, least-squares fit of those data for which 6* at 
X ? 0.2L was greater than 1.50k yields the correlation: 

5 = 0.8500 + 0.0219 (23) 

\ ' X r 0.2L 

Using this correlation, 5 » 0.883 when the theoretical value of the dis- 
placement thickness at x * 0.2L is equal to 1.50k. 

The ratio of the laminar displacement thickness to the tile misalign- 
ment height is presented in Fig. 35 as a function of x/L. Note that 



over the range of test Reynolds number at this surface temperature. Thus, 
it should not be surprising that, if the critical value of the displacement 
thickness at x = O.IL is 0.75k, then the critical value of 6* at x : 0.2L 


is 1.50k. 



CONCLUDING REMARKS 


The theoretical heat-transfer distributions are compared with 
experimental heat-transfer distributions obtained in Tunnel B at AEDC 
using a 0.0175 scale model of the Space Shuttle Orbiter Configuration 
for which the first 80% of the windward surface was rcughened by a sim- 
ulated tile misalignment. Data were obtained for a Mach number of 8 over 
a Reynolds number range (based on model length) from 1.862 x 10® to 
7.091 X 10® with surface temperatures from 0.114 T^ to 0.435 T^. For 
the geometries and for the flow conditions of the present wind-tunnel 
test programs, the following conclusions are made. 

1. A study of the surface-pressure distributions, the Mach number at the 
edge of the boundary layer, and the heat-transfer distributions indicates 
that calculations based on the Var Ent flow model are in good agreement 
with the laminar data for x >^0.2L. For x £ O.IL, the calculations based 
on the Mod Newt NSE flow model are in reasonable agreement with the data. 

2. The theoretical solutions indicate that thinning the boundary layer 
by surface cooling increased the nondimensionalized value of the local 
heat-transfer coefficient. The theoretical increase was due to the in- 
creased gradients and viscous dissipation which result when the boundary- 
layer thickness decreased. However, the coaxial surface thermocouples 
used to measure the heating rates were not calibrated to show the varia- 
tions of the magnitude predicted by theory. 

3. Tile misalignment did not significantly affect the heat-transfer 
rate in regions where the boundary layer was either laminar or turbulent. 

4. The heat-transfer distributions were used to determine the transition 
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locations over a range of Reynolds number, surface cooling, and 
tile-misalignment heiyiit. The transition location moved significantly 
forward as surface cooling at the higher Reynolds number thinned the 
boundary layer. Furthermore, relatively large tile-misalignment heights 
were required to cause significant movement in the transition location. 

The tile-induced flow perturbations promoted transition only when the 
theoretical value of 6* at x = 0.1 L was less than 0.75 k. This critical 
value of 6*/k was independent of the flow-field model used to calculate 
6* and was independent of the angle-of-attack for these wind-tunnel tests . 
When 6* at X = 0.1 L was equal to, or greater than, 0.75 k, the experiment- 
ally-determined transition locations for a tile-roughened, cooled model 
were within 20% of the reference, smooth-body transition location at the 


same free- stream conditions. 
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Figure 1. - Sketch of Space Shuttle Orbiter configuration 
(Rockwell Drawing VL70-000139) used in Tunnel 
B at AEDC. 
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Fiqure 2 . - Phntoqraph of model showinq randomly niisaliqned tiles. 








Figure 4. - Streamwise variation of the Mach number at the edge of the boundary layer. 
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Theoretical BLIMP distribution (ref. 19) x * Q.501L 

(both theoretical distributions for Re ■ 1.9 x 10°, T^ = 0.423 T^) 
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Figure 5. - Variation of the Mach number across the laminar boundary 
layer for x » 0.5L, o = 30®. 
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Figure 7. - Concluded. 
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Figure 10. - The effect of the surface temperature distribution on the 

heat transfer distribution. Re , • 1.9 x lO®, « = 30®. 
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Figure 11. - Concluded. 










Figure 13. - Concluded. 
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Figure 14. - The effect of surface temperature on the theoretical M^.ch 
number distribution in a laminar boundary- layer a": a = 30 
for Re , * 7.1 X 106 at X = 0.255L. 
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Figure 15. ~ Shddowgrdphs 
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Figure 16. - Continued. 
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effect of tile misalignment height on the heat-transfer 
tribution for o * 40*. 
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(a) a = 30“ 

Figure 21. - Transition locations as a function of the Reynolds number 
behind a normal shock wave, data from Tunnel B. 




















A fillec symbols are for smooth model 


A open symbols are for grit-roughened model 



Figure 26. - Transition locations as a function of the Reynolds num! 

behind a normal shock wave, data from Tunnel F, 3 = 30 
0.16 Ty, < 0-27 T^ (data from refs. 27 and 28). 
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(i) a = 30“ 

Figure 27. - The theoretical valup* of (ROg/M^) using the Var Ent flow 
model. " tr 
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Figure 28. - The relative transition location c, cc a function of 6 for 

Mod Newt NSE Theory, a = 30°. ^ 





Figure 28. - Concluded. 























Figure 32. - Concluded 









Figu 'e 33 • - Coi “luded. 














